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ABSTRACT 

We extend the MOND analysis to a sample of 17 high surface brightness, early-type 
disc galaxies with rotation curves derived from a combination of 21cm Hi line obser- 
vations and optical spectroscopic data. A number of these galaxies have asymptotic 
rotation velocities between 250 and 350 km/s making them among the most massive 
systems (in terms of baryonic mass) considered in the context of MOND. We find that 
the general MOND prediction for such galaxies - a rotation curve which gradually 
declines to the asymptotic value - is confirmed, and in most cases the MOND rotation 
curve, determined from the mean radial light and gas distribution, agrees in detail with 
the observed rotation curve. In the few cases where MOND appears not to work well, 
the discrepancies can generally be understood in terms of various observational errors 
- such as incorrect orientation angles and/or distances - or of unmodelled physical 
effects - such as non-circular motions. 

The implied mass-to-light ratios for the stellar disc and bulge constrain the MOND 
interpolating function; the form recently suggested by Zhao & Famaey (2005) yields 
more sensible values than the one traditionally used in MOND determinations of 
galaxy rotation curves. 

Key words: gravitation - dark matter - galaxies: structure - galaxies: kinematics 
and dynamics. 



1 INTRODUCTION 

As an acceleration-based modification of Newtonian dynam- 
ics, MOND makes general predictions about the contrasting 
form of rotation curves in low surface bri ghtness (LSB) and 
high surface brightness (HSB) galaxies (|Milgroml 1 19831 ) . If 
surface density is roughly proportional to surface bright- 
ness, then LSB galaxies have low internal accelerations and 
MOND predicts a large discrepancy between the visible and 
Newtonian dynamical mass within the optical image; in par- 
ticular, in these systems the rotation curves should slowly 
rise to the final asymptotic circular velocity. On the other 
hand, HSB galaxies, within the bright inner regions, are in 
the high acceleration, or Newtonian, regime. MOND pre- 
dicts that there should be a small discrepancy between the 
observable mass and the dynamical mass within the bright 
optical image; moreover, the rotation curves should rapidly 
rise and then fall in an almost Keplerian fashion to the final 
asymptotic value. Hi observations of several galaxies span- 
ning a range of surface brightness certainly exhibit this gen- 
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eral trend (|Casertano fc van" Gorkomlll99lh but the sample 
is small and the data are inhomogeneous. 

With respect to LSB galaxies the MOND pre- 
dictions have been spectacularly confirmed by Hi 
observations of suc h obj ec ts over the past decad e 
i|McGaugh fe de Blokl 1 19981: ISanders fc Verheiienl Il998l ). 
Most recently, Milgrom fc Sanders! ( 20061 ) have extended 
this analysis to several galaxies with extremely low mass or 
rotation velocity - systems where the observed distribution 
of neutral gas dominates the detectable mass budget of the 
galaxies and the mass-to-light ratio of the stellar component 
effectively disappears as a free fitting parameter. They 
demonstrated that MOND generally accounts for the form 
and the amplitude of rotation curves in these systems. 

At the other extreme, there has not been much con- 
sideration of the high luminosity or HSB galaxies. The 
gravitational force distribution probed b y planetary nebulae 
kinem atics in several elliptical ga laxies (iRomanowskv et al.l 
l2003h is consistent with MOND |Milgrom fc Sanderslbooj ) 
- primarily the basic prediction of a small discrepancy in 
high surface brightness systems. However, because of the 
uncertain degree of anisotropy in the velocity distribution 
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of these probes, rotation curves remain a far more unam- 
biguous tracer of the radial dependence of force. 

A sample of eight rapidly ro tating, but late-type, 
disc g alaxies has been presented by ISpekkens fc Giovanellil 
l|2006h . but these objects are apparently not as centrally con- 
centrated or bulge-dominated as the early type rapid rota- 
tors that we require in order to confront the generic MOND 
predictions. The rotation curves of several early type disc 
galaxies in the published lite rature at the tim e were consid- 
ered in terms of MOND by ISandersl l|l996h . but the data 
were quite inhomogeneous and the derived rotation curves 
were not of the same quality as in more gas rich systems. The 
basic problem is that high luminosity, HSB galaxies tend to 
be of early-type (SO-Sb) and have generally low gas content 
with patchy distributions. 

Recently, a more homogeneous sample of early-type 
disc galaxies with high-quality Hi ob servations has be- 
come available dNoordermeer et~aL I l2005l. hereafter N05). In 
iNoordermeer et aL (2007, hereafter N07), these data were 
combined with high-resolution optical spectroscopic obser- 
vations to derive rotation curves for 19 systems. These 
galaxies have morphological types ranging from SO to Sab, 
and have, for the most part, an apparent spheroidal bulge 
component and a high central surface brightness. They 
are mostly massive systems, with maximum rotation ve- 
locities of 200 - 500 km/s. Significantly, surface photom- 
etry in the B- and R-bands exists for all these galaxies 
(INoordermeer fc van der Hulstl2007l ); this is necessary to es- 
timate the distribution of visible matter and the resulting 
Newtonian acceleration. 

In this paper, we extend the MOND analysis to this 
sample of HSB, rapidly rotating galaxies. The rotation 
curves follow the general MOND predictions. In the bright 
inner regions there is little discrepancy between the observed 
and Newtonian rotation curves, while at larger radii the ro- 
tation curves often decline to an asymptotically constant 
value. Applying the MOND formula to the mass distribution 
as traced by the visible light and Hi emission, we find that 
in most cases the observed rotation curve is well matched by 
the MOND prediction. There are a few exceptions where the 
observed rotation curves appear to deviate from the MOND 
predictions, and these cases are considered in more detail. 

Because the accelerations range smoothly from the 
Newtonian to the MOND regime, this sample of rotation 
curves can be used to constrain the form of the MOND in- 
terpolating function \i. Using the typically applied form for 
this function we find that the mass-to-light ratios of the 
bulge and disc are often implausible; in particular, in sev- 
eral cases a higher M/L is required for the disc than for the 
bulge. Usi ng the simpler interpola ting function recently sug- 
gested by IZhao fc Famaevl (120061 ). the M/L values become 
more reasonable with the disc values generally comparable 
to or lower than those of the bulge. 

The remainder of this paper is structured as follows. In 
section O we describe the sample of galaxies and the obser- 
vational data. Section [3] describes the procedures we used to 
derive the predicted MOND rotation curves. In section [4] 
we present the general results, while individual galaxies are 
discussed in section \5\ Finally, in section [6] we draw some 
conclusions. 



THE SAMPLE AND OBSERVATIONAL 
DATA 



The rotation curves in lN07l are based on the Hi data from 
InoEI. which in turn were obtained in the framework of the 
Westerbo rk survey of Hi in spiral and irregular galaxie s 
(WHISP; iKamphuis et ailll996l : Ivan der Hulst et alj[200ll) 
A full description of the sample selection is given in iNOTl ; 
here we briefly mention a few points relevant for this study. 

All galaxies in have well-resolved Hi velocity fields 
(at least 5 to 10 independent beams across), and sufficient 
Hi across the disc to define a complete two-dimensional field. 
Moreover, their velocity fields show no obvious evidence for 
large-scale deviations from circular motions about the cen- 
tre of the galaxy: no strongly interacting or strongly barred 
galaxies were selected to avoid complications introduced by 
non-circular motions. 

Of the 19 galaxies in N07, two have been dropped 
here. The first (UGC 624) has an asymmetric velocity field, 
undermining the quality of its rotation curve. A second 
(UGC 4605) is nearly edge-on and no photometry is avail- 
able to trace the stellar mass distribution. 

The observational properties of the remaining sample 
of 17 galaxies is given in Table 1. Note that all distances are 
estimated from the systemic velocity, corrected for Virgo- 
centric inflow and assuming h=0.75. 

For all of these objects, the Hi velocity field has been 
supplemented by higher resolution single slit optical spec- 
tra in the central regions. This is evident as a much higher 
density of data points defining the inner observed rotation 
curves. While higher resolution data is certainly desirable in 
the region where the rotation velocity is changing rapidly, 
a note of caution is necessary. These data points are based 
upon single slit spectra and the implicit assumption is that 
the observed velocities along the slit can directly be trans- 
lated into rotational velocities (using the position angle and 
inclination as defined by the Hi velocity field) . But streaming 
motions of gas, for example, by weak non-axisymmetric dis- 
tortions in the central regions, can introduce uncertainties in 
the interpretation of optical velocities as circular velocities. 
Moreover, in any blind fitting procedure, the more numer- 
ous optical data points should not be given weight equal the 
outer Hi points, to avoid that the fitted model is determined 
entirely by the optical data. 

In addition to these considerations, there are the usual 
uncertainties which complicate the interpretation of the ro- 
tation curves as tracers of the radial force distribution. There 
are, in a number of cases, at least weak non-axisymmetric 
distortions (e.g. bars) in the inner region that introduce er- 
rors in the derived rotation curve. For other objects, the 
gas layer in the outer regions is clearly warped. In principle, 
these effects were taken into account as well as possible by 
using a tilted-ring model in the deriva tion of the rotation 
curves l|Begemanlll987l ; lBegemanlfl989h : a set of rings each 
with its own inclination, position angle and rotation velocity 
was fit to the radial velocity data (see |N07| for details) . But 
this procedure contains implicit physical assumptions - such 
as constant azimuthal velocity in a ring - which may or may 
not be realistic. The overall consequence of these effects is 
that the presented rotation curves cannot with certainty be 
taken as a true indicator of the force. The indicated error 
bars are an attempt to quantify these uncertainties, but we 
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Table 1. Sample galaxies: basic data. (1) UGC number; (2) alternative name; (3) morphological type; (4) distance 
(based on Hubble flow, corrected for Virgo-centric inflow and assuming h=0.75); (5) and (6) B- and R-band lumi- 
nosities (corrected for Galactic foreground extinction); (7) central surface brightness (R-band; corrected for Galactic 
foreground extinction); (8) maximum and (9) asymptotic rotation ve locity and (10) cen t ripeta l acceleration at last 
point in rotation curve. Column (3 ) was taken from NED, (4) from Noordcrm eer et al.l {2005) and (5) - (7) from 
Noordcrme er fc van der Hulstl |2007t) . 



UGC 
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(1) 
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(3) 
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(5) 
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(7) 


(8) 


(9) 


(10) 


2487 


NGC 1167 


SA0- 


67.4 


8.79 


10.47 


17.00 


390 


320 


0.46 


2916 




Sab 


63.5 


4.09 


3.37 


16.97 


220 


200 


0.29 


2953 


IC 356 


SA(s)ab pec 


15.1 


4.79 


5.50 


15.72 


315 


260 


0.44 


3205 




Sab 


48.7 


3.53 


2.99 


17.17 


240 


210 


0.40 


3546 


NGC 2273 


SB(r)a 


27.3 


1.59 


1.84 


15.84 


260 


190 


0.43 


3580 




SA(s)a pec: 


19.2 


0.328 


0.311 


17.87 


127 


120 


0.20 


3993 




SO? 


61.9 


1.85 


1.84 


16.99 


300 


250 


0.41 


4458 


NGC 2599 


SAa 


64.2 


5.55 


5.86 


16.00 


490 


260 


0.34 


5253 


NGC 2985 


(R')SA(rs)ab 


21.1 


3.44 


3.05 


15.56 


255 


220 


0.31 


6786 


NGC 3900 


SA(r)0+ 


25.9 


1.47+ 


1.50 


15.88 


230 


210 


0.48 


6787 


NGC 3898 


SA(s)ab 


18.9 


1.56 


1.72 


15.35 


270 


250 


0.64 


8699 


NGC 5289 


(R)SABab: 


36.7 


0.964 


1.05 


16.47 


205 


180 


0.45 
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NGC 5533 


SA(rs)ab 


54.3 


4.79 


5.92 


16.27 


300 


230 


0.17 


11670 


NGC 7013 


SA(r)0/a 


12.7 


0.745 


0.879 


15.59 


190 


160 


0.40 


11852 




SBa? 


80.0 


2.33 


2.17 


17.50 


220 


160 


0.09 


11914 


NGC 7217 


(R)SA(r)ab 


14.9 


2.00 


1.84 


16.05 


305 


290 


3.48 


12043 


NGC 7286 


SO/a 


15.4 


0.160 


0.107 


18.71 


93 


85 


0.18 



+ No B-band data available in Noordcrmeer & van der Hulst (2007); Lb taken from LEDA. 



cannot exclude the possibility that in individual cases, the 
true circular velocity deviates from the observed rotation 
curves. 

Finally, we use the pho tometric data from 
iNoordermeer fc van der Hulsti (|2007h to derive the contribu- 
tion to the gravitational field from the stars. Most of these 
galaxies show evidence for two distinct components in the 
light distribution - a central spheroidal bulge and a more 
extended disc component. INoordermeer fc van der Hulstl 
presented bulge-disc decompositions, where a bulge with a 
generalised Sersic profile was subtracted from the total light 
distribution to give the remainder as the disc component. 
Assuming that the disc is a flat system, with a vertical scale 
height one fifth of the radial scale length, and comparing the 
ellipticity of the isophotes in the central bulge-dominated 
regions with those of the outer disc, they could also esti- 
mate the intrinsic axial ratio of the bulges. Note that the 
bulge-disc decomposition introduces additional uncertainty 
in our mass models. Usually, a range of bulge and disc 
parameters can be derived to fit the observed intensity 
distribution. We have verified that, as long as the bulge 
and disc mass-to-light ratios are comparable, differences 
in the bulge-disc decompositions have little effect on the 
mass models. However, when the mass-to-light ratios differ 
strongly, this is no longer the case, and errors in the 
decomposition propagate into the mass models. 



3 PREDICTED MOND ROTATION CURVES 

The procedure here for generating a MOND rotation 
curve is basically the sam e as described previously (e.g. 
ISanders fc McGaughlliooi l. One complication for this sam- 
ple is the presence of a bulge which may have a different 
mass-to- light ratio than the disc component. This adds un- 



certainty due to the decomposition procedure as well as a 
second free parameter. 

Basically one takes the photometry in each of the two 
components, bulge and disc, as a faithful tracer of the stel- 
lar mass distribution (mass-to-light ratio in each component 
does not vary with radius) . Then given the estimated intrin- 
sic shape of the bulge and assuming that the disc is a flat 
system, one may calculate, using the usual Poisson equa- 
tion, the radial distribution of force (assuming a M/L ratio 
fo r each component). This procedure is described in detail 
in lNoordermeerl l|2006l ). The observed surface density of neu- 
tral hydrogen is included (multiplying by a factor of 1.3 to 
account for primordial helium) assuming this component to 
be in a disc with thickness equal to that of the stellar disc. 

The Newtonian force qn is then converted into a "true 
force" , g by using the MOND formula 

gfj,{g/a ) = g N . (1) 

Here ao is the critical acceleration (the one ne w parameter of 
the t heory), which we fix at 1 x 10 _8 cm/s 2 (|Bottema et al.l 
2002). is the MOND interpolating function which must 
have the asymptotic dependence fj,(x) — -> x when i < 1 (the 
MOND limit) and jj,(x) — > 1 when x 3> 1 (the Newtonian 
limit). Unlike the case of the LSB galaxies, where |<?|/ao < 1 
everywhere, in these objects there is a gradual transition 
from the Newtonian limit in the inner regions to the MOND 
limit in the outer parts. Therefore, one might expect the 
form of /x to play some role in the fitted models, at least 
in the derived M/L values. Two different forms have been 
considered in the recent literature. Most often one takes 

Ms) = -7=5, (2) 
VI + x 

but more recently IZhao fc Famaevl (|2006l ) have suggested a 
simpler form which is consistent with a multi-field relativis- 
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tic extension of MOND l|BekensteiiJl2004 ): 



iFamaev et al.l (|2006T > have demonstrated that both functions 
produce reasonable fits to observed rotation curves, but the 
second (eq. 3) generally yields lower mass-to-light ratios of 
the visible components. Here we will apply both forms and 
consider the plausibility of the inferred M/L ratios as a con- 
straint on fj,. 

The circular velocity is determined from the true force 
(v 2 /r — g) and the M/L ratios of the two components are 
adjusted to achieve the optimal fit to the observed curve. 
Thus, for the objects in this sample, where there is a clear 
bulge component, there are two free parameters in the fitting 
procedure. In view of the uncertainties associated with the 
errors on the observed rotation curves (see previous section) , 
a simple x 2 -niinimisation does not always result in the best 
possible fit. Instead, we have used the x 2 -statistics only as 
a starting point for inspection by eye of the individual fits, 
and in some cases adjusted the mass-to-light ratios by hand 
to achieve a better fit to the overall shape of the rotation 
curve. 

Finally, it should be noted that all dis tances were taken 
to be the Hubble law values given by |N05l (assuming h=0.75 
and corrected for Virgo-centric inflow). Because most of 
these galaxies are quite distant, Hubble expansion is gener- 
ally a fair indicator, and distance was not allowed to be an 
additional free parameter in these fits. However, a few galax- 
ies lie more nearby, in which case the Hubble distance can 
be off by a considerable amount. In such cases, the quality 
of the fit might be improved further by allowing the distance 
to vary within reasonable limits (see also section [5| . 



4 RESULTS 

The resulting MOND rotation curves are compared to the 
observed rotation curves in figure [1] Here the points with 
error bars show the observed velocities. The thin, black 
solid line indicates the Newtonian rotation curve of the 
bulge component, while the dashed and dotted curves show 
the contributions of the stellar and gaseous discs respec- 
tively. These Newto nian curves are, in form, those given by 
iNoordermeerl ([2006) . The bold, blue solid curve is the re- 
sulting MOND rotation curve assuming the interpolating 
function given by eq. 3. The implied mass-to- light ratios of 
the two stellar components are given in Table 2. 

The error bars shown here require a word of explana- 
tion and a note of caution. There are three sources for the 
indicated errors: 

First are the formal velocity fitting errors. For the ve- 
locities from the optical spectra, these simply come from 
the derivation of the centre of the emission line. For the ve- 
locities from the Hi data, they are given by the non-linear 
least square routine that fi ts tilted rings to the entire two - 
dimensional velocity fields (Bcgcman 1987; Bcgcman 1989). 
In the tilted ring method it is assumed that deviations from 
planar circular motion can be modelled as pure rotation in a 
ring with a different inclination or position angle on the sky. 
While this may be a fair assumption in the outer regions 
where warps are often seen, it is most likely not correct in 



the inner regions where systematic deviations in the veloc- 
ity field may result from non-circular motion due to a non- 
axisymmetric gravitational field. In such cases, the formal 
errors underestimate the true uncertainties. 

Secondly, the tilted ring analysis was performed sepa- 
rately for the approaching and receding sides of the galaxy 
and this often resulted in different rotation curves for the 
two different sides. This is also an indication of asymme- 
tries in the velocity fields, but not the sort of bi-symmetric 
asymmetries that would be expected from bar-like distor- 
tions. The error on the rotation curve due to these effects is 
taken to be one-fourth the difference in the rotation curves 
of the two sides. 

Thirdly, there is an error resulting from the uncertainty 
in the inclination. This was estimated by eye looking at the 
scatter in the fitted inclination resulting from the tilted ring 
program, and is most often the largest contribution to the 
error budget. 

The first two sources of errors were added quadratically 
to produce the error bars shown in figure [1] The latter un- 
certainty is indicated by the gray shaded bands. Combined, 
they should be taken as a rough indication of the uncertainty 
on the derived rotation curve, but, again, they are not true 
statistical error bars. This means that statistical tests, such 
as evaluating different models by comparing \ 2 ar e meaning- 
less and misleading. One should also recall that the points 
on the rotation curve are not all independent since the Hi 
velocity field was sampled every half beam-width of the Hi 
observations. 

Figure [1] shows that many of the observed rotation 
curves in this sample of galaxies exhibit a rapid rise fol- 
lowed by a gradual decline to a constant asymptotic value 
(as in e.g. UGCs 2487, 2953, 3993, 4458). This is a funda- 
mental prediction of MOND for galaxies with high central 
surface brightness systems. Moreover, for these extreme cen- 
trally concentrated galaxies, the Newtonian rotation curves 
generally coincide with the MOND curve for the inner few 
kpc (i.e. the mass discrepancy in the inner regions is small) . 
This is particularly evident in UGC 11914 where the entire 
measured rotation curve is in the high acceleration regime, 
and, as MOND would predict, the discrepancy between the 
Newtonian and the observed curve is small. 

Considering this sample of rotation curves overall we 
see that in about 10 out of the 17 cases the MOND rota- 
tion curve closely agrees with the observed rotation curve. In 
some cases (e.g. UGC 2953) the MOND curve is a spectac- 
ular reproduction of the observed curve. In five other cases 
(UGCs 3205, 3580, 6786, 9133 and 11670) MOND correctly 
predicts the general behaviour of the observed curve but 
misses details. For two galaxies (UGC 6787, UGC 11852), 
there are significant differences between the MOND predic- 
tion and the observed rotation curve, in the overall shape as 
well as the details. This is comparable to previous samples 
(see Sanders & McGaugh 2002) where, in roughly 10% of 
the cases a substantial mismatch between the MOND pre- 
dictions and the observed rotation curves was found. 

Given the uncertainties in the derivation of rotation 
curves and their precision as a tracer of the radial force 
distribution, such an occasional mismatch is not surprising 
and does not pose an insurmountable problem for MOND. 
Nonetheless, we discuss these individual problematic cases 
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Figure 1. MOND fits to the observed rotation curves of early-type disc galaxies. Data points show the observed rotation curves and 
error bars give the combined uncertainties due to the measurement errors and kinematical asymmetries. The grey shaded bands give the 
allowed range due to inclination uncertainties. Thin black solid, dashed and dotted lines give the contributions from stellar bulge, disc 
and gas respectively. The thin green (grey) line gives the Newtonian sum of the individual components and the bold blue (grey) lines 
gives the total MOND rotation curve, assuming the interpolating function fi from eq. 3. 



and the possible sources of the disagreement in the following 
section. 

With respect to the implied mass-to-light ratios, the 
standard interpolating function (eq. 2) often requires values 
which are implausibly high, in particular for the disc (e.g. 
UGCs 3993, 6786). This problem is particularly apparent 
when the disc M/L-ratios are compared to those of the bulge: 
in many cases the former is larger than the latter (e.g. 2953, 



3993, 5253, 6786), whereas naively, one would expect the 
redder bulge to have a higher M/L. With the simpler form 
given by eq. 3 however, the fitted mass-to-light ratios of the 
disc are lower, and usually comparable to or smaller than 
the bulge M/L. This indeed suggests that the more gradual 
transit ion between New tonian and MOND dynamics implied 
by the IZhao fc Famaevl interpolating function is superior. 

In figure[2] we compare the fitted global (i.e. bulge+disc) 
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Figure 1 - continued 



mass-to- light ratios (with p given by eq. 3) t o the predictions 
of rec ent population synthesis models of iBell fc de Jond 
(|200lf ) and lPortinari et al.1 ( j2004T >. The values generally scat- 
ter above the models, but if we ignore the two highest points 
(UGC 3993 and UGC 11914), the overall trend is similar 
- the redder the galaxy, the higher the M/L. It is also of 
interest that these two discrepant objects have the largest 
uncertainties in inclination. There is a systematic offset to- 
ward higher M /L values than would be predicted by the 
IBell fc de Jond mo dels, but, as we see in comparison with 
the Portinari et al.l models, this can be accommodated by 
adjustments in the initial stellar mass function. 



5 INDIVIDUAL OBJECTS 

UGC 3205: The MOND rotation curve deviates significantly 
from the observed curve in the inner 10 kpc. This may be a 
result of streaming motions in the weak bar in this galaxy, 
which cause the observed velocities to deviate significantly 
from the local circular velocity and lead to asymmetries 
in the observed Hi velocity field |N05). Outside the region 
where the velocity field is influenced by the bar, the observed 
rotation velocities agree well with the MOND predictions. 

UGC 3580: The MOND rotation curve is generally sim- 
ilar to the observed curve but the match can be improved 
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Figure 1 - continued 



considerably by increasing the distance to the galaxy by 
about 20%. This is within the uncertainty in the distance 
to this relatively nearby object. 

UGC 3993: The overall form of the MOND rotation 
curve agrees reasonably well with the observed rotation 
curve; however, the M/L ratios for the bulge and disc seem 
unacceptably large with both the Zhao-Famaey interpolat- 
ing function (> 8 for both component s) and certainly with 
the standard function (26 for the disc) . iNoordermeerl (|2006T ) 
notes that this galaxy deviates strongly from the R-band 
Tully-Fisher relation in the sense of having a luminosity that 
is too low for its rotation velocity. The galaxy is nearly face- 



on (assumed inclination of 20 degrees) which means that 
there is a large uncertainty in the amplitude of the derived 
rotation curve. This is a plausible source of the anomaly: 
the galaxy is not under-luminous but the amplitude of the 
rotation curve may have been overestimated due to an incli- 
nation over-correction. In this case, the mass-to-light ratios 
of the bulge and disc will be correspondingly lower as well. 

UGC 6786: In the MOND rotation curve shown, the 
M/L values of the bulge and disc have been adjusted to 
match the decline in the observed rotation curve in the outer 
regions. In this case, the MOND rotation curve is then too 
high in the inner regions. There are large side-to-side differ- 



8 



R.H. Sanders & E. Noordermeer 



Table 2. Fitted mass-to-light ratios in the R-band for discs and bulges separately, 
and for the entire galaxies. Columns 2-4 refer to the fits with the 'standard' 
function fi(x) = x/Vl + x 2 , whereas columns 5-7 are for the form suggested by 
IZhao fc Famaevl l|200rj l: fj,(x) = x/(l + x). 
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Figure 2. Fitted global (bulge+disc) mass-to-light ratios in the 
R-band vs. B-R colour. Red (gray) and black data points show the 
mass-to-light ratios from th e fits with the standard interpolating 
function fi and the one from I Zhao fc Famaevl (2006) respectively. 
Errorbars show the photometric errors and the formal fitting er- 
ro rs on M/L. The bold d iagonal line shows the relation derived 
by Bell &: de Jon 3 ll200ll) , while the thin lines show models from 
iPortinari et ail (|2004|7 for different assumptions for the initial stel- 
lar mass function. The dashed regions show linear extrapolations 
beyond their bluest and reddest models. 



ences in the inner rotation curve - of about 100 km/s; such 
asymmetries complicate the interpretation of the observed 
radial velocity as a rotation velocity in the inner regions. 
However, at radii larger than 4 kpc the rotation curve is 
quite symmetric and the Hi and Ha data are consistent. 



This is an extreme bulge-dominated galaxy; it appears to 
be more like an elliptical with a small disc rather than a 
normal Sa. 

UGC 6787: Here the MOND rotation curve differs in 
form as well as in detail from the observed rotation curve. 
This galaxy also deviates strongly from the R-band Tully- 
Fisher relat ion in having too la rge a rotation velocity for its 
luminosity |Noordermeer| [2006h . A possible problem here is 
the distance estimate. This galaxy lies in the direction of 
the Ursa Major cluster but has a redshift 200 km/s larger 
than the high velocity envelope of this cluster. It is possible 
that the true distance is larger than assumed here, but that 
the galaxy has a lower apparent Hubble velocity because 
of in-fall into the UMa cluster. However, to achieve reason- 
able agreement with the MOND rotation curve, the distance 
must be increased by almost a factor of two, which seems 
implausible. Furthermore, changing the distance does not 
improve MOND's ability to reproduce the observed 'wiggles' 
in the rotation curve. An unmodelled warp beyond 16 kpc 
could cause the apparent rise in the rotation curve, but then 
the inclination of the gas layer in the outer regions would 
need to approach 90° (edge on) to produce the apparent rise 
of almost 30 km/s. The dip in the rotation curve between 
12 and 18 kpc coincides with a change in the fitted position 
angle of the tilted rings, indicating deviations from simple 
iso-planar gas motion. This does imply that some of the ap- 
parent structure in the rotation curve may be due to the 
tilted ring fitting algorithm. 

UGC 9133: Although the general form of the MOND 
rotation curve agrees with the observed curve, the detailed 
structure is not reproduced well. As in the case of UGC 6787, 
the fluctuations in the derived rotation curve of this galaxy 
are coincident with significant fluctuations in the fitted tilted 
ring parameters; in particular there is a shift of more than 
10° in the position angle of the inclined rings coincident 
with the bump in the rotation curve at about 60 kpc. Here 
the tilted ring program is trying to accommodate a real sys- 
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tematic twist in the two-dimensional velocity field - a twist 
which may well be due to non-circular motion, in addition to 
variations in the orientation of the gas disc. In any case, ar- 
tifacts in the rotation curve which are too dependent upon 
the tilted ring fitting program should be interpreted with 
caution. 

UGC 11670: The inner dip in the rotation curve (at 
about 4 kpc), which is not reproduced by MOND, coincides 
with a fairly dramatic shift in the fitted position angle of the 
tilted rings. The true culprit is likely to be non-circular mo- 
tion in the potential of the bar. The velocity field is rather 
messy in the outer regions; thus, the observational uncer- 
tainties on the outer three points on the rotation curve may 
be larger than indicated with the error bars. The MOND 
fit to the rotation curve is considerably improved if the dis- 
tance is increased by 30% - a possibility because this is the 
closest galaxy in the sample. 

UGC 11852: This galaxy is clearly warped. There is a 
large shift in the fitted inclination of the tilted rings coincid- 
ing with the apparent decline in the rotation curve. Again, 
it is hard to exclude that possibility that the warp is accom- 
panied by non-circular motions. In fact, IN07I reported large 
residual velocities with respect to the tilted ring model, in- 
dicating that some streaming motions are indeed present. 
A reasonable fit to this galaxy's rotation curve can also be 
achieved if it lies at about 75% of its adopted Hubble law 
distance. 



6 CONCLUSIONS 

It is rare when an astrophysical hypothesis makes actual 
predictions in the realm of extragalactic phenomenology and 
rarer still when those predictions are verified. But with re- 
spect to the systematic form of galaxy rotation curves for 
low surface brightness galaxies, the modified Newtonian dy- 
namics has been exceptional in this respect. I Milgrom] pre- 
dicted in 1983, before a substantial population of low LSB 
galaxies was discovered, that in such objects the rotation 
curve should rise throughout the visible disc to the final 
asymptotic value. Since that time, a population of these ob- 
jects has been discovered and observed, and this behaviour 
is confirmed. 

A supplementary prediction concerned high surface 
brightness galaxies: in these objects the rotation curve 
should, after a rapid rise in the very central regions, slowly 
fall and approach the constant asymptotic value. 

The sample presented here offers the chance to test this 
prediction. These early-type galaxies are mostly of high cen- 
tral surface brightness with a substantial bulge component. 
The rotation velocities are among the largest observed, rang- 
ing up to 500 km/s. Such galaxies have been previously 
difficult to observe because of a relatively low and patchy 
gas content, but with the upgraded Westerbork array and 
supplementary optical spectroscopic observations in the in- 
ner regions, it has become possible to map two dimensional 
velocity fields and derive rotation curves that extend well 
beyond the bright inner regions. The general prediction of 
MOND with respect to such high surface brightness objects 
is confirmed. 

MOND is now seen to successfully account for the sys- 
tematics and details of galaxy kinematics from very low ro- 



tation, or acceleration, systems l|Milgrom fc Sanders! lioolih 
to the galaxies presented here with extremely high rotation 
velocities, all with the same value of the critical accelera- 
tion parameter ao- This extends the range of viability of 
MOND over a range of a factor of 10 in rotation velocity 
and 10 4 in luminosity /mass. It therefore becomes increas- 
ingly problematic to interpret the success of MOND as it 
merely being an algorithm for explaining the relation be- 
tween dark and visible matter, as some have suggested. A 
more natural explanation is to see MOND as a manifesta- 
tion of a fundamental modification of gravity or inertia in 
the low acceleration regime - that MOND is the key to new 
gravitational physics. 
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